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Abstract 

Peripheral heavy ion reactions at ultra relativistic energies have large an¬ 
gular momentum that can be studied via two particle correlations using the 
Differential Hanbury Brown and Twiss method. We analyze the possibili¬ 
ties and sensitivity of the method in a rotating system. We also study an 
expanding solution of the fluid dynamical model of heavy ion reactions. 


1 Introduction 

Collective flow is one of the most dominant observable features in heavy ion re¬ 
actions up to the highest available energies, and its global symmetries as well as 
its fluctuations are extensively studied. Especially at the highest energies for pe¬ 
ripheral reaction the angular momentum of the initial state is substantial, which 
leads to observable rotation according to fluid dynamical estimates [T|. Further¬ 
more the low viscosity quark-gluon fluid may lead to to initial turbulent instabilities, 
like the Kelvin Helmholtz Instability (KHI), according to numerical fluid dynamical 
estimates (2], which is also confirmed in a simplified analytic model [3]- These 
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turbulent phenomena further increase the rotation of the system, which also leads 
to a large vorticity and circulation of the participant zone one order of magnitude 
larger than from random fluctuations in the transverse plane |4][5l[6j. 

The Differential Hanbury Brown and Twiss (DHBT) method has been introduced 
in 13 . The method has been applied to a high resolution Particle in Cell Relativistic 
(PICR) fluid dynamical model |10| . 


2 The two particle correlation 


The pion correlation function is defined as the inclusive two-particle distribution 
divided by the product of the inclusive one-particle distributions, such that m- 
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where pi and p 2 are the 4-momenta of the pions and k and q are the average and 
relative momentum respectively. 

We use a method for moving sources presented in Ref.[l4j. In the formulae the 
h = 1 convention is used and k and q are considered as the wavenumber vectors. 
The correlation function is: 


C{k,q) = l + 
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where 
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Here R{k,q) can be calculated jl4] via the function and we obtain the R{k,q) 
function as 

Rik,q) = Re[J{k,q)J{k,-q)] (4) 

The corresponding J{k,q) function will become 
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For the phase space distribution we frequently use the Juttner (relativistic Boltz¬ 
mann) distribution, in terms of the local invariant scalar particle density the Juttner 
distribution is (15] 
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where C„ = Airm^TK 2 {m/T). We assume a spatial distribution: 


G{x) = ‘jn{x) = 'jUs exp 
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Here rig is the average density of the Gaussian source, s, (or fluid cell) of mean 
radius R. 

Asymmetric Sources: we have seen in few source model examples [9] that 
a highly symmetric source may result in correlation functions that are sensitive to 
rotation, however, these results were not sensitive to the direction of the rotation, 
which seems to be unrealistic. We saw that this result is a consequence of the 
assumption that both of the members of a symmetric pair contribute equally to the 
correlation function even if one is at the side of the system facing the detector and 
the other is on the opposite side. The expansion velocities are also opposite at the 
opposite sides. The dense and hot nuclear matter or the Quark-gluon Plasma are 
strongly interacting, and for the most of the observed particle types the detection 
of a particle from the side of the system, - which is not facing the detector but 
points to the opposite direction, - is significantly less probable. The reason is partly 
in the diverging velocities during the expansion and partly to the lower emission 
probability from earlier (deeper) layers of the source from the external edge of the 
timelike (or spacelike) FO layer. 

For the study of realistic systems where the emission is dominated by the side 
of the system, which is facing the detector, we cannot use the assumption of the 
symmetry among pairs or groups of the sources from opposite sides of the system. 
Even if the FO layer has a time-like normal direction, the {k^a^) factor yields a 
substantial emission difference between the opposite sides of the system. 

The correlation function, C{k,q) is always measured in a given direction of the 
detector, k. Obviously only those particles can reach the detector, which satisfy 

k^&fj, > 0. Thus in the calculation of C{k,q) (see Fig. 1) for a given k- direction 
we can exclude the parts of the freeze out layer where k^a^ < 0 (see Eq. (10) 
of Ref.jlb] or Ref.jl 8 |. For time-like FO a simplest approximation for the emission 
possibility is Pesc(a^) cic k^^u^{x) m- 


3 The DHBT method and fluid dynamical re¬ 
sults 

Based on the few source model results the Differential HBT method [9) was in¬ 
troduced by evaluating the difference of two correlation functions measured at two 
symmetric angles, forward and backward shifted in the reaction plane in the partic¬ 
ipant c.m. frame by the same angle, i.e. at 77 = ±const., so that 


AC{k,q) = C{k+,qout) - C{k_,qout)- ( 8 ) 

For the exactly zLx -symmetric spatial configurations (i.e. k+x = and 
k+z = —k-z), e.g. central collisions or spherical expansion, AC{k, q) would vanishi 
It would become finite if the rotation introduces an asymmetry. 

The sensitivity of the standard correlation function on the fluid cell velocities 
decreases with decreasing distances among the cells. So, with a large number of 
densely placed fluid cells where all fluid cells contribute equally to the correlation 
function, the sensitivity on the flow velocity becomes negligibly weak. 
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Figure 1: (color online) The dependence of the standard correlation function in 
the direction from the collective flow, at the final time. From ref. pH] . 


Thus, the emission probability from different ST regions of the system is essential 
in the evaluation. This emission asymmetry due to the local flow velocity occurs 
also when the FO surface or layer is isochronous or if it happens at constant proper 
time. 



q (1/fm) 

Figure 2: (color online) The differential correlation function AC{k^q) at the 
final time with and without rotation. From ref. m- 

We studied the fluid dynamical patterns of the calculations published in Ref. 
[21, where the appearance of the KHI is discussed under different conditions. We 
chose the configuration, where both the rotation [I], and the KHI occurred, at 
b — O.lbmax with high cell resolution and low numerical viscosity at LHC energies, 
where the angular momentum is large, L « lO^h [13]. Fig. |^shows the DHBT for 
the FD model. 
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The standard correlation function is both influenced by the ST shape of the 
emitting source as well as its velocity distribution. The correlation function becomes 
narrower in q with increasing time primarily due to the rapid expansion of the system. 
At the initial configuration the increase of |fc| leads to a small increase of the width 
of the correlation function. 

Nevertheless, in theoretical models we can switch off the rotation component 
of the flow, and analyse how the rotation influences the correlation function and 
especially the DCF, AC{k,q). 



q (1/fm) 

Figure 3: (color online) The Differential Correlation Function (DCF) at average 
pion wavenumber, k = 5/fm and fluid dynamical evolution time, t = 3.56fm/c, 
as a function of the functions of momentum difference in the ’’out” direction 
q (in units of 1/fm). The DCF is evaluated in a frame rotated in the reaction 
plane, in the c.m. system by angle a. From ref. [10) . 

Fig. [^compares the standard correlation functions with and without the rotation 
component of the flow at the final time moment. Here we see that the rotation 
leads to a small increase of the width in q for the distribution at high values of |fc|, 
while at low momentum there is no visible difference. 

In F\g^AC{k,q) is shown for the configuration with and without rotation. For 
k — 5/fm the rotation increases both the amplitude and the width of AC. The 
dependence on |fc| is especially large at the final time. 

Fig. shows the result where the rotation component of the velocity field is 
removed. The DCF shows a minimum in its integrated value over q, for a = 
— 11 degrees. The shape of the DCF changes characteristically with the angle a. 
Unfortunately this is not possible experimentally, so the direction of the symmetry 
axes should be found with other methods, like global flow analysis and/or azimuthal 
HBT analysis. 

Finally we separated the effect of the rotation by finding the symmetry angle 
where the rotation-less configuration yields vanishing or minimal DCF for a given 
transverse momentum k. This could be done in the theoretical model. We did this 
for two different energies, Pb+Pb / Au+Au at ^/snn = 2.36/0.2 TeV respectively. 
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Figure 4: (color online) The differential correlation function (DCF) with and 
without rotation in the reference frames, deflected by the angle a, where the 
rotation-less DCF is vanishing or minimal. In this frame the DCF of the original, 
rotating configuration indicates the effect of the rotation only. The amplitude 
of the DCF of the original rotating configuration doubles for the higher energy 
(higher angular momentum) collision. From ref. [10] . 


while all other parameters of the collision were the same. The deflection angle of the 
symmetry axis was a = —11/ — 8 degree^ respectively. In these deflected frames 
we evaluated the DCF for the original, rotating configurations, which are shown in 
Fig.[^ This provides an excellent measure of the rotation. 


4 Summary 

We show that two particle correlation measurements can be sensitive to the rotation 
of the emitting system. The analysed model calculations show that the Differential 
HBT analysis can give a good quantitative measure of the rotation in the reaction 
plane of a heavy ion collision. 
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